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1. Introduction 

Despite tremendous progress in our understanding of fundamental strings in the past 
decade, we are still very far from a single quantitative prediction to be observed in exper- 
iments! The main reason for this unsatisfactory state of affairs is that the natural scale 
in which the stringy effects become important (Planck scale) is much smaller than the 
scale we can probe in high energy scattering experiments. However, it has been known for 
a while that cosmology provides a window to test fundamental physics, since the scales 
of particle physics theories become relevant as the universe evolves in its early history. 
It is our belief that the most likely area for a confrontation between string theory and 
experiments is through extracting cosmological consequences of string theory which may 
be measurable today. Many of the riddles of astrophysics, for example 'dark matter' and 
large scale structure of the universe may turn out to be related to fundamental strings. It 
is with such a program in mind that we study string cosmology in this paper. 

To see the effects of strings we have to trace back the universe to its earliest times. 
Here already string theory provides us with a surprise: strings in very small spaces behave 
the same way as strings in very large spaces [1,2]. This is true at least in the simplest 
known cases (toroidal backgrounds) and is suspected to be a general principle of string 
theory, because strings should be 'blind' to scales smaller than their own natural scale. 
Thus the universe does not start at a zero size, as that would be equivalent to infinitely 
large universe, rather we should start with a universe of the size of the string scale, i.e., the 
Planck length. This a — > 1/a duality (in Planck units) already hints of a major modification 
to Einstein equations at early times, as ordinary theory of general relativity is not invariant 
under such a transformation. It has been suggested that, roughly speaking, string theory 
necessitates the existence of two kinds of matter fields, "dual" to each other. For large radii, 
one type becomes massless and dominates, while for small radii the other type becomes 
massless and dominates. This suggests an effective "doubling" of the number of space 
dimensions (see e.g. [2,3]). Needless to say to have an exact treatment of dynamics in 
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such a situation is beyond our scope at the present time. This was one of the main reasons 
why the cosmological ideas in [2] were rather difficult to test, as the dynamical equations 
governing string cosmology were not properly understood. 

It turns out that it is possible to introduce a consistent dynamics which describes 
slowly evolving states of strings in a duality invariant manner. This involves, in particular, 
taking into account the dilaton field which transforms under the duality [4-7]. 

It is a fundamental consequence of string theory that the gravitational interaction 
is effectively described not by the metric alone but by the coupled system of the metric 
and the dilaton. Dilaton of course should get a mass to avoid a scalar component of 
gravitational interaction of massive bodies but this probably happens at a later stage of 
evolution of the Universe. That implies that one cannot not ignore the dilaton dynamics 
at early "stringy" phases of cosmological evolution. As we shall see, the assumption that 
dilaton is constant is in general inconsistent with the cosmological equations describing 
evolution of matter with the "stringy" equation of stated. 

The main body of this paper is devoted to the discussion of the duality invariant cos- 
mological equations describing the evolution of the "scales" of the metric and the dilaton. 
We shall study their solutions emphasizing departures from Einstein's gravity theory. For 
example, we will find that a constant energy density not only does not lead to inflation, 
but it rather slows and eventually stops expansion!!. Similarly, we will discover that if 

1 Dilaton dynamics was ignored in most previous discussions of string cosmology at finite 
temperature (see e.g. [8] and refs. there). Dilaton was not set equal to constant in ref. [9] where, 
however, the space was not taken to be compact and hence many of the stringy effects, and in 
particular duality invariance was not discussed. 

2 We shall describe the expansion/contraction in terms of the original unrescaled ( "<r -model" 
) metric. Though the effective gravitational constant is then time dependent , this is the metric 
string directly interacts with and hence the one measured by stringy "rods". Similar point of 
view was expressed in [10]. The use of this metric is appropriate at early phase of cosmological 
evolution when dilaton is still massless (for a discussion of a possible inflation in string cosmology 
described in terms of the rescaled metric after the dilaton already got mass see [11]). Let us note 
also that the use of the "unrescaled" metric is particularly natural in the duality symmetric case 
since it has a simple transformation law under the duality. 
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the winding modes (by "winding modes" we mean a more general notion of stringy states 
whose masses grow with a growing scale of the universe) are not annihilated (which may 
happen if they are not in thermal equilibrium) they halt the expansion. Using our equa- 
tions we can study some simple aspects of the cosmological scenario proposed in [2] in 
which the three dimensionality of extended space was suggested to be a consequence of 
string cosmology. 

The equations we consider apply to both critical and non-critical strings and we study 
them both in the presence and in the absence of matter. The simplest solution correspond- 
ing to non-critical strings was first proposed by Myers [12] as a way to change the dimension 
of target space in string theory and was then used for construction of cosmological solutions 
in [13]. These solutions were further extended by Mueller [14] who found non-trivial cos- 
mological solutions for which the radii of the toroidal universe change with time. We will 
show that specializing his solution to d=2 one recovers (a region of) the two dimensional 
black-hole solution, which has been found to be related (in the leading order approxima- 
tion) to a conformal theory [15]!. We find that the more natural setup to discuss d=2 
cosmology is to go to an infinite fold cover of SX(2, R)/U(l). This will avoid getting time- 
like loops and give us infinitely many universes parametrized by an integer n, such that 
the future of particles in the universe n could be in any universe m > n. The exact metric 
suggested for the d=2 Euclidean black hole [17] which leads to a puzzle within a black-hole 
interpretation seems perfectly consistent once we use such a cosmological interpretation 
(corresponding to the universal cover of SL(2, R)/U(l)). 

The organization of this paper is as follows. In section 2 we derive the basic equations 
of string cosmology in the adiabatic approximation and discuss their duality invariance 
property. In section 3 we specialize to critical strings and consider some basic elements of 
string cosmology by studying the solutions of these equations. We find a simple intuitive 

3 As was already pointed out in [5], the euclidean continuation of the d=2 Myers solution 
coincides with the euclidean solution independently discovered in [16]. 
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way to think about the solutions based on the analogy with a mechanical model of a particle 
rolling on a potential under the influence of a damping force. In section 4 we concentrate 
on some aspects of non-critical string cosmology and its relation to the d—2 black hole. We 
also briefly discuss how string one loop effects might modify the non-critical cosmological 
solutions. In section 5 we present our conclusions. In appendix A some aspects of the 
solutions to our cosmological equations are worked out. 

2. General Structure of Cosmological Equations 

Let us consider strings propagating in a large universe of spatial dimension N. If N is 
not the critical spatial dimension of strings (i.e., if N < 25 (9) for (super) strings), we could 
have one of two possibilities: either we have some N c of compact internal dimensions which 
make up the rest of the spatial dimensions which we take to be static and to correspond 
to conformal theories on the worldsheet, or we are dealing with a "non-critical" string. In 
other words, if we set 

c=^(d c -N-N c ) 

where d c is the critical spatial dimension, in the first case c = 0, and in the second case 
c 7^ 0. For simplicity we take the iV— dimensional space to be a periodic box (torus) of 
lengths a,i = expAi ( i runs from 1 to N). However, many of the ideas developed in this 
paper can easily be extended to more general cases. 

It is a well known principle in string theory that as long as we concentrate on slowly 
varying fields, we can hope to describe the global aspects of the dynamics of the the- 
ory reasonably well by concentrating on massless modes, and keeping only the leading 
derivative terms ('adiabatic' approximation). The types of massless modes depend on a 
particular choice of vacuum backgrounds strings are propagating in, and on whether we 
are considering bosonic strings, superstrings, or heterotic strings. But they always include 
a gravitational field and a dilaton (for simplicity we shall ignore the antisymmetric tensor 
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field). If the radii are large compared to the string scale the leading order terms in the low 
energy expansion of the tree level gravitational-dilaton effective action of the closed string 
theory are [18] (we shall absorb the gravitational coupling constant into <p an d set o/=l)§ 



Note that e^ plays the role of coupling 'constant' in string theory. The first question we 
would like to address is whether we should trust this action at all for Planckian size spatial 
dimensions, i.e., when ~ 0. The answer is no, if we wish to take all the modes in the 
above action seriously. The above action manifestly breaks the duality symmetry of string 
theory. The reason is that when we approach Planckian size universes the momentum 
modes (which correspond to inhomogeneity in the metric and dilaton fields) become mas- 
sive, and the winding modes which wrap around the box and were ignored previously will 
have to be taken into account, as they have mass comparable to that of the momentum 
modes. This 'doubling of the space' seems very difficult to deal with, and on top of that 
we have to remember that the oscillator (higher level) modes are also going to become 
important since they will be of the same energy as the momentum modes. 

However, the following observation allows us to still take part of the above action se- 
riously. The point is that the space independent parts of the above fields have zero energy 
independent of the size of the universe, and thus we can naturally concentrate on them, 
setting the other Planckian modes to zero. Note in particular that these spatially indepen- 
dent modes are common to both the momentum and winding excitations, so truncating to 
these modes alone will be consistent with duality. Therefore we should expect that even 
though the above action naively seems to pick momentum modes over winding modes, 
i.e., the space over the dual space, concentrating on spatially constant fields restores the 
duality symmetry. Thus we should trust the above action as far as the overall scales of 

4 The volume of the internal space, if there is any, can be gotten rid of by a shift in <j>. 




(2-1) 
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the universe are concerned and not so much as a description of inhomogeneities of the 
universe. And in addition we should assume that the fields are varying slowly with time. 

Let us see how the duality invariance <Zj — > l/a^ is restored in the above action once 
we concentrate on zero modes. Let us therefore consider the metric and the dilaton field 
of the form 

N 

ds 2 = -dt 2 + af(t)dx 2 , 
i=i 

a l = e x ^ , <P = <P(t) . 

It is useful to introduce the "shifted" (and rescaled by 2) dilaton field (f which absorbs the 
space volume factor 

N 

i=i 

Then the action becomes 



r , N 

S = - dt e-^y^o [c-G 00 £A? + G°V 

J i=l 

We have kept Gqq as we need to vary the action with respect to it in order to get the full 
set of equations of motion ( in the equations of motion G 00 is set to —1). Dot denotes the 
time derivative. The above action is manifestly invariant under the duality transformation 
[4-6] 

Xi -> -A; , 4>—>4>-Xi , if — > ip (2.2) 

(here % is any number among 1,..,N ; more general duality transformations are obtained 
by combining these basic ones). 

The above action is written in the absence of stringy "matter" . To discuss cosmology 
it is important to incorporate matter into the system. Assuming that the effective string 
coupling is small let us take the matter to be a gas of (free) string modes in a thermal 
equilibrium at the temperature Then matter contribution to the action is represented 
by 

S m = J dt^/^GVo F(Xi, Pv^GVo) ■ 



Here F is the (one loop) free energy and can be represented in terms of the one loop string 
partition function in a torus of radii e Xi and periodic Euclidean time of perimeter (3^— Goo 
, Z = —f3F. The adiabaticity assumption implies that we can replace constant radii and 
(3 by functions of time. Taking the full action 

r - 

S = S Q + S m = - dt y^G^ [ e -^( c - G 00 J2^ + G 00 ^ 2 )- F(Xi, /3a/^Goo)] 

"* i=l 

and varying it with respect to \i,<f> and Goo and we find the following equations^ 

N 

c-J2^i+^ = ^E , (2.3) 

i=l 

\ i -<p\ i =\e*P i , (2.4) 

N 
i=l 



where 



p = = _ dF_ 

5Xi d\i ' 

E is the total energy of the matter and P$ is the pressure in the i-th direction times the 
volume. These equations imply a modified conservation law for the energy (following from 
the invariance under reparametrizations of time ) 

N 

1=1 

Since F = F(X(t), (3(t)) this is equivalent to the conservation of the entropy S = (3 2 dF/d(3, 
i.e. our matter is indeed evolving adiabatically. This means that for a given radius deter- 
mined by Aj the temperature f3~ x adjusts itself so that S remains constant. Solving the 



5 Similar equations were discussed in [9] and also in [6] where the matter was represented by the 
energy-momentum tensor of a macroscopic string. Even though we derived the above equations 
using the canonical ensemble, they are valid even when the canonical ensemble description breaks 
down (as is the case near Hagedorn temperature) and we have to use the more fundamental 
microcanonical ensemble. 
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adiabaticity condition we can express (3 in terms of A^ and hence represent E as a function 
of Xi alone, 

E(X) = E(X,(3(X)) . 

In what follows we shall consider only this new function E(X). Since entropy is not changing 
we find 

Note that the time derivative of (|2.3|) vanishes once we use (|2.4|) and ( |2.5|) , and so the 
non-trivial content of equation ( |2.^| ) is that it fixes the constant of integration in terms of 
c. The above equations are duality invariant since F is invariant under A — > — A for a fixed 
temperature (E(X) and /3(A) are also invariant and Pi changes sign under the duality). 
Rewritten in terms of the original dilaton field <fi equations (2.3)-(2.5) take the form 

N N 

c-£A? + (20-5>) 2 = e 2 *p , (2.7) 



A, - (20 - *i) Ai = k^Pi , {2.1 



i=i 

N N 

20-^A l -^A l 2 = ie^p (2.9) 
i=i j=i 

where p = E/V , pi = Pi/V and V = expJ^Ai is the space volume. To compare these 
equations with standard cosmological equations in which the dilaton is ignored consider 
setting cj) (which is related to the coupling constant) to be constant. Then the equation 
(|2.7|) takes the well known form (with flat space) , 

(;) , = * a = ^+ " < 2 - 10 > 

for some constant G, where —c/(N 2 — N) plays the role of the cosmological constantly. 
However, it is easy to see that if c = a solution with <fi = const is consistent with all the 



6 It is amusing to note that in this setup the problem of the vanishing cosmological constant 
seems to be related to the question of being on or off-criticality in string theory. This connection 
(which was already mentioned in the past) seems worth further study. 



three equations only if 

N 
i=l 

i.e., for a matter with the vanishing energy momentum trace, as is effectively the case 
for a gas of massless particles in a ther mo dynamical equilibrium. This "radiation -type" 
condition is not satisfied in a high temperature phase of string thermodynamics and is the 
reason for stringy departure from Einstein's gravity. 

Let us discuss some general aspects of equations ( p73| ) , (|2.4|) and (|2.5|) . It is important 



to notice that the equations for (|2.4|) are the same as that of a point particle in the 
presence of a time dependent potential \e v E(\) and damped or boosted (depending on 
whether (p is negative or positive) by a dilaton "friction" force. Note that within this 
mechanical interpretation the "integration constant" c in (2.3) plays the role of a fixed 
energy of the system. 

Also, the equation for (p can be summarized by subtracting equation (|2.5| ) from ( |2~3| ) 
and defining y = e _VP , which leads to 

y + cy= l -E . (2.11) 

This equation for y can be interpreted as that for the endpoint of a spring, with spring 
constant c, and with the external force ^E(X). These interpretations allow us to develop 
an intuitive picture for the nature of the solutions to our equations. As it turns out, the 
behavior of solutions will be rather different depending on whether or not we are describing 
critical strings, i.e., whether or not c = 0. 



3. Critical String Cosmology 

This case corresponds to setting c = 0. The important qualitative simplification 
that happens in this case is that since E is positive (assuming no Casimir - like negative 



contribution to energy) from equation (|2.3|) we conclude that (p can never become zero. In 
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other words, (p never changes sign, and thus it provides a damping or a boosting effect in 
equation (|2.4| ) for all time. We will consider the damping case <p < 0; the boosting case 
can be obtained by time reversal. A strong reason to consider (p < is that otherwise 
the boosting effect of dilaton on A will invalidate the adiabatic approximation used in the 
derivation of the above equations. Also, growing (p together with expanding A implies the 
growth of the effective coupling exp (p in contradiction with a weak coupling assumption. 

From equation ( |2.5| ) it follows that Cp > and thus <p is growing in time. Since (p 
is negative and it will never cross zero we conclude that (as long as E ^ 0) (p continues 
growing towards and approaches as t — > oo. 

For simplicity let us consider the isotropic case , i.e. assume that all Xi are equal to 
each other and denoted by A so that a = exp A is the cosmological scale. The equations 
we get can be written as 

c-NX 2 + ip 2 = e^E , (3.1) 
\-ip\= \^P , (3.2) 
Cp - NX 2 = \^E , (3.3) 

where P = —N~ 1 dE/d\ and for critical strings we simply set c = in ( |3.1| ). 

In order to solve this system of equations we will need to specify E(X) as well as 
provide the initial conditions for A, (p and (p. E(X) encodes string thermodynamics by 
specifying how the total energy in the box has to change as a function of A in order to keep 
the entropy constant. Motivated by string cosmology, the function E(X) was introduced 
and studied using string thermodynamics in [2] ( some properties of E(X) were further 
clarified in [19], see also [20] ). Its basic structure is indicated in fig. 1. 

This function is duality symmetric under A — > —A. Its structure is well understood 
near A ~ and for very large A. In the first region near A ~ 0, which we call the 
Hagedorn region, the temperature of the strings is very close to the Hagedorn temperature 
and the energy is almost independent of A. It is given by E = ThS where Tjj is the 
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Hagedorn temperature (fixed for a given string theory) and S is the total entropy. For large 
enough A temperature drops significantly compared to the Hagedorn temperature, and the 
massive modes of string go out of equilibrium so that we are left with the massless modes. 
In this 'radiation dominated' region we thus have the usual relation of how E depends 
on the radius of the universe for a gas of massless particles. In fact, the temperature 
drops as 1/a = exp(— A) and E oc T N+1 a N oc a -1 . In other words E = Ce~ x where 
C is determined by S and the number of massless modes. The behaviour of E in the 
intermediate region between the Hagedorn region and the radiation dominated region is 
only partially understood but it has been shown that E is monotonically decreasing (in 
the Hagedorn regime) with increasing |A| [19], and it is natural to believe that this is 
true for all A. So if we start with the initial conditions that cp < the decreasing of the 
potential E with increasing |A| will force A to slide down the potential towards increasing 
values of |A|, until we emerge at the radiation dominated era. The universe will continue 
its expansion, and in fact it will approach the radiation dominated era of the standard 
cosmology. In other words, independently of the details of E, or the initial conditions, one 
finds that 20 = <p + N\, which determines the string coupling constant, rapidly approaches 

2 

a constant and the radius ultimately grows as a ~ A t N + 1 for large t, as is the case for 
standard radiation dominated cosmology. This is discussed in appendix A. 

What determines N, the number of extended spatial dimensions? One idea to explain 
why N < 3 was suggested in [2] . It was based on the possibility that if the space expands 
in more than three dimensions, the assumption of maintaining thermal equilibrium may 
not be a good one, as for example the winding modes of strings will have a hard time 
finding and annihilating one another. Thus they would fall out of equilibrium and will be 
around. It was suggested there that if winding modes are around this may stop expansion. 
The universe would not expand until it learns the lesson that it is only possible to expand 
in three or smaller number of dimensions. 

Let us note that Einstein's theory of gravity leaves us uncomfortable with the above 
suggestion that winding modes may prevent expansion. Recall that cosmological expansion 
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rate for a flat universe with vanishing cosmological constant satisfies A 2 



Gp in FRW 



cosmologies. In particular any form of matter, since it contributes to p, helps accelerate 
expansion. This in particular implies that the winding modes accelerate expansion in 
standard gravity theory, rather than prevent expansion. On the other hand the idea in [2] 
was based on duality which states that what usual matter does for the size of the universe 
(i.e. accelerates expansion) the dual matter (the winding states) should do for the dual 
size of the universe (i.e., accelerate contraction). This was a puzzle which was not fully 
resolved in [2] . We will now see that our equations resolve this paradox by showing that 
they drastically differ from Einstein's equations in such cases and indeed are consistent 
with the idea that winding modes slow down and ultimately stop the expansion. 

If the winding modes are around the energy E is going to grow with A, as larger A 
mean bigger boxes and thus more stretched winding strings with higher energies. The mass 
of the winding states will increase with A as expA. E growing with A means that A < as 
A will have to climb the potential E. So this slows the expansion and in fact ultimately 
stops it as it is shown in appendix A. More generally, it is shown there that if we take for 
large A 



then for a positive the universe reaches a maximum size and then contracts back. Con- 
stant matter energy density corresponds to a = N and it prevents expansion even more 
strongly! So we conclude that dilaton suppresses the inflationary mechanism based on con- 



radiation dominated era where a = —1). 

The duality symmetry of E implies that while the winding modes (providing the 
growth of E at large A) oppose expansion, the momentum modes (providing the growth of 

7 As was already noted in the Introduction, this does not rule out the possibility of inflation 
at a later stage of evolution when dilaton gets a mass and it is appropriate to use the redefined 



E oc exp(aX) 



stant energy density^. If 




a is negative the universe expands forever (as is the case for the 



metric. 
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E at large negative A ) oppose contraction. As a result, in the absence of equilibrium the 
universe sees the effective E represented as a dashed curve in fig. 1 and the radius of the 
Universe will be oscillating between maximal and minimal values. 

So the basic picture one is led to is a universe which oscillates in many directions 
for a while around the Planck scale (maybe within a few orders of magnitude), until by 
coincidence it starts expanding in smaller number of directions (the most likely one being 
the largest possible dimension consistent with maintaining equilibrium, namely 3). Then 
it expands forever and we find ourselves in the radiation dominated era of the standard 
cosmology. This may also explain the large entropy problem. The entropy may be large 
because during this "trial and error" period of the early universe, as the strings were out 
of equilibrium we would generate ( if this process takes long ) a lot of entropy. In other 
words, the large entropy problem may be related to having an 'old' universe. 

An important issue to resolve in string cosmology is to explain the absence of massless 
dilatons at the present time. One should expect that somehow a potential is generated for 
the dilaton and it will have to sit at its minimum picking up a mass. The mechanism of 
how this precisely should happen in string theory is not known yet and is the biggest gap 
in connecting string cosmology to observable cosmology. This is presumably related to the 
fundamental question of how supersymmetry breaking takes place in string theory while 
maintaining vanishing cosmological constant. 



4. Non-critical String Cosmology 

The case of non-critical string cosmology is an analog of non-vanishing cosmological 
constant case of cosmology in standard theory of gravity. We shall consider the case of 
c > 0. Cosmological solutions for the toroidal space and non-zero c in the absence of 
matter ( E = ) has been previously studied in [14]. One obtains from (|2.11|) ( setting 
E = ) 

y = e~ v = A sin 2bt , b 2 = c/4 . 
13 



From (|2.4j ) ( with Pj = 0) one learns that 



Xi = ki , ki = const 
which can be easily integrated to yield 

Aj = Xio + qi In tan bt . 
subject to the condition resulting from equation (|2.3|) 

N 

E«. 2 = 1 • 



i=l 



These solutions can also be extended to c < [14] . Note that for c > these solutions 
are singular, i.e., the fields blow up at finite time. In fact, as it is easy to see from ( |2.11|) , 
the oscillatory nature of solution forces y, which should be positive, to cross zero in finite 
time ( so that ip blows up in finite time) . This could in principle be balanced by a positive 
energy density appearing on the right hand side of (|2.11[) . 



For the above singular solutions the adiabatic approximation which went into their 
derivation breaks down, and so we should not trust them near the singularity. Even so, 
one would expect that there should be a solution for any reasonable initial condition, if 
we believe in completeness of string theory. This leads us to expect that there should be 
an exact conformal theory which asymptotically reproduces Mueller's solutions. For the 
case iV = 1 we will now see directly that this is indeed the case. In fact the corresponding 
conformal theory turns out to be the SL(2, R) /U(l) coset model which was recently linked 
to the black hole geometry in [15]. 

For N = 1 we have qi = ±1. Let us take the plus sign for q as the other sign (giving 
the dual solution [5]) corresponds simply to shifting of time. We thus have 

tp = cpo — In sin 2bt , A = Ao + hi tan bt , b 2 = 4 . 
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This background describes a universe with one spatial dimension with the metric given by 
(we fix a particular value for Ao and assume that x is periodic with period r , e.g. 2tv ) 

ds 2 = -dt 2 + e 2X dx 2 = -dt 2 + b~ 2 tsm 2 bt dx 2 . 

The universe starts at t = with zero size and grows to infinite size at t = n/2b. Let us 
change our coordinates. Let 

u = sin bt e x , v = sin bt e~ x . (4.1) 

We find that 

,9 , _o dudv 

ds 2 = -b 2 , 

1 — uv 

and that the unshifted dilaton field <fi is given by 

2(p = tp + A = 2</> - 21n cos bt = 2</> - ln(l - uv) . 

This is the same metric and dilaton background found in [15] in the region < uv < 1 (here 
uv = sin 2 6t), i.e. the region between the horizon and singularity. There is one difference 
though: here we are working on a periodic space, so x is identified with x + r. This means 
that we have to identify 

(w, v ) ~ (cm, a~ l v ) (4-2) 

i.e. we get a wedge in the region between horizon and singularity (see fig. 2). Of course 
the system of equation we looked at is also valid for infinite radius r, which implies that we 
can in fact reproduce the full region between horizon and singularity with no identification. 
In this case in order to complete the metric we will have to add regions / and JiT as well, 
and we end up getting the full black hole solution. It is amusing that the transformation 
( fl.2|) is an exact symmetry of the conformal theory (in an appropriate basis it corresponds 
to conjugating SL{2,R) group element by the diagonal matrix (a 1 / 2 , a -1 / 2 )). So we can 
consider an infinite orbifold of SL(2, R)/U(l) by the group generated by this symmetry 
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and obtain the compact universe solution. One should note that since u = v = is a fixed 
point of this transformation we end up getting a new singularity at u = v = 0, which in 
our cosmological interpretation corresponds to the initial time t = 0, where the universe 
has zero radius. There is another wedge branch of the solution in regions I and III which 
has periodic time which touches the cosmological region at the singular point u = v = 0. 

Note that in this language the duality of conformal theory [17, 21] corresponds (for 
appropriate choice of radius) to the standard a — > 1/a duality of our one dimensional 
universe, i.e. bt — > \ — bt changes uv — > 1 — uv and A — > —A. 

Another aspect of Mueller's solution is that we can continue time past the singularity 
at t = tv /2b. In particular, if we introduce 

u\ = —sin bt e x , v\ = —sin bt e~ x 

we get an identical copy of the interior region of black hole (on the overlap u\ = —u , v\ = 
—v) where the universe now starts at the lower singularity, and evolves till t = ir/b at which 
it reaches zero size, it grows as t increases until it becomes infinite in size again at t = 3iv/2b. 
Then again we use the coordinate transformation 

U2 = sin bt e x , V2 = sin bt e~ x 

where the radius shrinks to zero size at t = 2n /b. By now the universe has undergone two 
oscillations, but we have covered SL(2, R) only once. For example, if we take x = 0, the 
corresponding element of SL(2, R) , as we evolve in time, is represented by the matrix 

/ cos bt sin bt \ 
y —sin bt cos bt J 

(it is known that the SL{2, R)/U(l) has two copies of the interior region of black hole-to 
get only one we could have used PSL(2, R)/U(l)). But now we will not identify (1*2,^2) 
with (u, v), as that would have given us periodic time. Instead we continue forever. This 
means that in the above matrix realization of SL(2, R) we are not identifying t — > t + 2n /b. 
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In other words the cosmological interpretation suggests that we go to an infinite fold cover 
of SL(2, R). In this way we end up with a universe which had no beginning and no end, 
and it undergoes infinitely many oscillations. Again in this picture we have the option of 
choosing a finite radius or infinite radius for our space. If we have a finite radius, then 
we will have to restrict to the wedge region between the horizon and singularity shown in 
fig. 2, otherwise we will have to add the regions I and III to complete the metric. In this 
latter interpretation an observer in region / or III of one of these universes, may decide 
to enter the next universe. The Penrose diagram for this series of universes is shown in 
fig. 3. It is amusing to note that in the Euclidean version of each of these universes we have 
a periodic coordinate which suggests a thermal bath interpretation. This suggests, even 
in the cosmological picture, we get radiation from the horizons. It would be interesting to 
study this further. 

Now we ask the question if this one dimensional cosmological solution should be taken 
seriously. Naively the answer would be no, because we have violated the assumption of adi- 
abaticity of the fields near the singularity, when the size of the universe is infinitely large. 
However, since we know there is an exact conformal theory with the correct asymptotic 
behavior, we are led to expect that there exist a corrected metric which gives the exact 
answer. Based on comparison with the coset model an exact metric for the Euclidean black 
hole was suggested in [17]. As was checked in [22] the corresponding metric - dilaton back- 
ground solves the a model conformal invariance conditions in the 3-loop approximation. 
When one analytically continues the conjectured metric to the region between horizon and 
singularity it reads 

ds 2 = -dt 2 + b~ 2 — & * dx 2 , q = 2/k = 8/9 . 

l + qtan 2 bt ' 1 

If this metric is indeed exact, it suggests that the space-time is actually nonsingular at 
t = tx /2b so that we see no singularity in the metric! In the finite radius scenario (where 
we introduce an additional singularity at u = v = 0) we have a universe starting at zero 
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radius, reaching a finite maximum radius at finite time and then contracting back to zero 
radius and continuing this indefinitely. The conjectured form for dilaton tp is unchanged 
from the leading approximation and is thus still singular at t = n/2b. The cosmological 
interpretation of the black hole gets rid of a puzzling consequence of this conjectured exact 
metric: if analytically continued to the region behind black hole singularity the metric 
becomes Euclidean between 1 < uv < 9 with a singularity at uv = 9. This is avoided in 
our interpretation because the region beyond the singularities are never reached (they are 
connected to universes labeled by ±oo). 

It would be interesting to find exact conformal field theory generalisations of higher 
dimensional Mueller-type cosmological solutions!. The coset theories like SO(d : 1) / SO(d— 
1, 1) which should be analogs of the (anti) DeSitter backgrounds (see Bars [21] and refs. 
there) does not, however, correspond to a flat iV-space. The generalization of the leading 
order conformal invariance equations to the case of isotropic homogeneous metric with a 
curved space (with positive, zero or negative curvature k ) is given by (cf. (3.1)-(3.3) ) 

c - NX 2 + p 2 = -kN(N - l)e" 2A ' 

X - (pX = -k(N-l)e~ 2X ' 
(p - NX 2 = 

(the absence of the correction in the third equation is due to the fact that the "potential" 
in the present case is "classical", i.e. does not depend on the dilaton). Let us assume that 
c < 0, k > 0. Then at large negative A (small times) the solution is approximately given 
by the DeSitter metric and a constant dilaton <fi = ^((p + NX). In this limit, however, the 

8 In such a case it would be natural to allow the full moduli of toroidal compactifications to 
become involved. The picture is a more or less straight-forward generalization of the equations 
we have considered, and we would simply have to study geodesic equations on the corresponding 
moduli spaces of toroidal compactifications. 
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adiabatic approximation is not fully reliable. At large positive A the effect of the space 
curvature becomes irrelevant and the solution approaches the Mueller's isotropic solution 

A = Ao + iV _ 2l n tanh bt — > const , <p — > — bt . 

So the universe is born att = Oasa<i = iV + l DeSitter space and evolves into the product 
of a time line and an iV-sphere at large t. This appears to be a direct higher dimensional 
(Minkowski signature) analog of the d = 2 (euclidean) "cigar" metric. 

So far we have discussed non-critical string cosmology in absence of matter. It is 
natural to ask how can one introduce 'thermal' non-critical matter, in order to have a 
richer non-critical string cosmology. It turns out that this can be done at least formally. 
The usual way we consider thermal ensemble in field theory is by making the time Euclidean 
with period f3. For non-critical strings the Liouville field plays the role of time and this 
suggests making it Euclidean and periodic with period (3. The computation of thermal 
free energy, which is a one loop computation generalizing that of [23] can be carried out 
explicitly for the non-critical strings c m < 1. For c m = 1 noncritical strings at radius a, 
for example, one finds that the necessary modular integral has been evaluated in [24] in a 
completely different context. This is related to the S 1 x S 1 partition function integrated 
over the fundamental domain and turns out to be the same as that for the critical N = 2 
strings partition function [25] ( or, equivalently, for c m = 1 non-critical N = 2 string 
partition function). One finds for the free energy 

(3F = ln[r 2V 2 (r)r 1 2 (f)] + (r^p) 

where r = ir 2 = i(3a/87r, p = i(3/8ira and rj is the Dedekind eta function. Solving the 
adiabaticity condition one can obtain the energy E(a) from the above expression for F. 
For large entropy and small a we have 

E(a) = 1 



a + a 1 
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where A depends on the entropy. The profile for E(a) is very similar to that of critical 
strings (fig.l) except that there is no flat Hagedorn region and we only get a maximum 
at A = 0. This is to be expected since at c m = 1 there is no exponential degeneracy of 
states, as there is only one massless scalar (and some additional discrete states [26]) and 
so there is no limiting temperature. In fact it has been suggested [27] that introduction 
of temperature for non-critical strings might lead to a more clear phase diagram picture 
for non-critical strings. Note that for large a, E(a) behaves like the energy for that of a 
massless particle (E oc 1/a). For yet larger radii, the temperature drops to zero, but E 
does not go to zero. It becomes negative and behaves as 

This is a Casimir - type contribution which is independent of the temperature. It is in 
accord with the result of [23] which corresponds to computing the free energy at zero 
temperature. Note that this Casimir effect is absent in critical (super)-strings, as the 
cosmological constant vanishes at one loop. Putting the potential E(a) on the right hand 
side of our equations will modify Muellers solution. The intuition based on the A rolling 
on the potential hill E(a) implies that the solutions are still singular (and blow up in finite 
time) . 

These ideas can be applied to c m < 1 as well. In this case we have a one dimensional 
field theory (corresponding to the Liouville field), which we can identify with a space 
coordinate x. The relevant equation to solve for the dilaton at the tree level is c — <p' 2 = 
which gives rise to the familiar linear dependence of the dilaton: ip = —\fc x. However, 
the one loop partition function [23] modifies the above equation to 

c-ip' 2 = -he v 

where h is a positive constant depending on which minimal model we are dealing with. 
Solving the above equation we find 

i ■ 1,2 Vex 

ip = ipo — In sinh . 

2 

This suggests a truncation of the Liouville space to a half-line. Again we are going be- 
yond the validity of adiabatic approximations but we are getting the hint that the linear 
dependence of dilaton is modified at one loop. 
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5. Conclusions 

We have discussed some aspects of string cosmology taking into account the dynamics 
of the dilaton field which is needed to restore the duality symmetry in string theory. We 
have considered both critical and non-critical strings, the non-critical strings being the 
analog of having a non- vanishing cosmological constant in the Einstein theory. 

For critical strings we find a major modification of the description of the early universe. 
In particular, winding string states, if not annihilated, can halt the expansion of the 
universe. This is against the intuition based on standard gravitational theory where a 
constant energy density leads to inflation. We have seen that string dynamics is consistent 
with the following picture of cosmology which needs to be verified. A nine dimensional 
spatial universe of Planck size which we take to be a periodic "box" , at the Hagedorn 
temperature, expands in all directions. Not being able to get rid of winding modes the 
expansion stops and the universe contracts to even smaller than the Planck scale. Now the 
momentum modes cannot annihilate one another and the contraction stops. The universe 
starts to grow once again. The universe oscillates for a while in this fashion. In this way 
one can in principle generate a lot of entropy. Fluctuations cause different directions to 
expand at different times. One of these fluctuations leads to three dimensions expanding, 
and now the winding modes can get annihilated and the universe expands, getting rid of 
the winding modes. As the temperature drops, the Planckian string states get suppressed 
by the Boltzman factor and we are left with the massless modes in a three dimensional 
expanding universe, which looks very much like the radiation dominated era of the standard 
cosmology. It is remarkable that the radiation dominated era of the standard cosmology 
emerges with no fine tuning. 

How can one test this picture? One idea is to treat strings classically as one does for 
cosmic strings. So we can study the evolution of a universe filled with these strings which 
interact with each other by cutting and rejoining. The main difference with the cosmic 
strings is that we have a different evolution equation for the universe, and the coupling 
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constant for the cutting and rejoining which is related to dilaton field is also evolving. This 
approach will presumably answer questions of entropy production and homogeneity of the 
early universe in the cosmological scenario described above. 

We have also studied non-critical string cosmologies, which seem qualitatively rather 
different from the critical string case. In particular we have related the empty universe 
cosmology in d = 2 to the two dimensional 'black hole' solution (more precisely to its 
universal cover). Could it be that stringy black-holes in 4d behave in a similar fashion? In 
particular, as we cross the horizon we enter another universe instead of encountering the 
singularity? It would be remarkable if the "cosmology-black hole" connection discovered 
in two dimensions has higher dimensional analogs. Indeed, the fact that the singularity of 
standard 4d black hole solution is also spacelike, i.e. appears at a given time rather than 
space, suggests a cosmological interpretation. 

To discuss non-critical string cosmologies with matter we have to introduce the notion 
of temperature for non-critical strings. We have made a suggestion of how this may be 
done (by a formal wick rotation of Liouville field and making it periodic with period f3 
as one does for c m = 25 in order to define critical string thermodynamics )i. Using this 
definition we have computed the thermal properties of d = 1 matter (which can be easily 
generalized to d < 1 as in [23]). It is an interesting question to see if there exists a matrix 
model analog of this finite temperature interpretation (for example are unitary matrix 
models related to thermalizations of standard matrix theories?) . 

We would like to thank M. Bershadsky, A. Guth, S. Jain, J. Polchinski, M. Tsypin, 
G. Veneziano and E. Verlinde for interesting discussions. A.T. would like to acknowledge 
J. Bagger for hospitality at Johns Hopkins University and various kinds of help. The re- 
search of A.T. was supported by NSF grant PHY-90-96198 and that of C.V. was supported 
in part by Packard Foundation and NSF grants PHY-89-57162 and PHY-87-14654. 

9 At higher genus due to the coupling of Liouville field to the curvature, this may be consistent 
only if P is quantized. 
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Appendix A. 

In this appendix we discuss the solution to the equations (|3.1|),(|3~2l),(|3~^) with c = 

and 

E oc exp(a\) (A.l) 

starting with initial conditions tp < 0, A > 0. In particular we show that if the total energy 
grows with volume, i.e., if a > the universe stops expanding in finite time and starts 
contracting, i.e., we reach A = in finite time as A is decelerating. If a < the universe 
expands forever and if a = the universe comes to a halt, without turning back. We will 
show in addition that if a = —1, as would be the case for radiation dominated era, the 
radius eventually grows as in standard Einstein's gravity in the radiation dominated era. 

One can reach the above conclusions more or less immediately by using the ideas 
developed in section 3. As discussed there starting with tp < 0, we conclude that tp cannot 
become zero, and thus it will always remain negative. Since tp > from ( |3.3|) we conclude 
that tp approaches zero and gets there at t = oo. Equation ( |3.2| ) implies that we can treat A 
as the position of a particle rolling on a potential described by e ip E(X)/2N with a damping 
term (because of tp terms in ( ft.2| )). The potential is getting weaker as its overall size is 
modulated by e^ which is decreasing because tp < and the damping term gets weaker 
as tp is approaching zero. For a < it is clear that since E is decreasing with A, A will 
continue to increase indefinitely. For a > 0, since the energy increases with A, the universe 
cannot grow indefinitely, and will reach a maximum and start contracting. This is almost 
obvious, were it not for the fact that the energy is modulated with e v which is getting 
smaller. So the universe could conceivably continue expanding at ever slower rates. This 
we will show is not the case. We shall find by a careful analysis of these equations that 
the universe reaches a maximum size and starts contracting after that. The fact that for 
a = (i.e., constant energy), the universe reaches a maximum size and stops is also clear 
from (|3.2j ) as the potential is flat in this case, and the damping term will eventually stop 
A. 
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First we note that P = -N~ l dE/d\ = —aE/N. Then we can eliminate e^E from 
equations (|3.1|) , (|372l) , ( |373| ) and obtain 



Cp-N\ 2 = -(-NX 2 + ip 2 ) . 
Now we can represent these equations in the first order form by defining 

I = A , f = <p 

to obtain 

/ = — + lf-^ 
• Nl 2 f 2 

The initial conditions are / < and / > 0. Actually we have another restriction, which 



comes from positivity of E in ( |3.1| ) 



\f\ > \VNl 



The important point to note in solving ([A.2 ) is that these equations are homogeneous. In 



other words, if we rescale (I, f) — > r(Z, /) and rescale t —> r 1 t we get a new solution. This 
means, in particular, that if we look for solutions with 

we get straight lines in the (/, I) plane which pass through (0, 0). Studying these particular 
solutions will give us a handle to study a qualitative behavior of all solutions. These straight 
line solutions can be easily found by noting that df jdl = f/l. Solving ( |A.3| ) we find three 
solutions 

I a 
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By the rescaling argument, / and / can approach zero only along the above lines. So in 
the region a > 0, if we start with the initial conditions Z>0, / < , |/| > VNl , the 
only possibility for not crossing the I = line is for the flow to approach zero through the 
line / jl = —\/N. But a simple check of the equation shows that this line is repulsive, i.e., 
the flows are driven away from it. They will be attracted to the next fixed line which is 
in the / < region. They will cross the / = line which means that at finite time A = 0. 
According to (|37|) A is not zero so the universe bounces back and A starts to decrease. 

For a = 0, the solution with I = is an attractive solution and all solutions will 
approach it. This simply means that if we start with A > the expansion will eventually 
stop because of the damping term. In fact, it is easy to write the general solution in this 
case: using ( |2.11|) and (|3.2|) ( with P = ) we get 

_ v E t 2 NA 2 1 t - 2VnA/E 

e v = — — , A = An H ^ln = . 

4 E y/N t + 2^/NA/E 

This solution approximately describes the evolution in the Hagedorn region (in fact, E is 
not quite constant and decreases with increasing A so the expansion never quite stops). 

If — y/~N < a < 0, the line / // = N/a is an attractive solution and all solutions will 
approach it. Note that for a = — 1 which corresponds to the case of the radiation dominated 
era, the fixed line is f/l = —N. For the fixed line solution 2</> = ip + NX = f + N I = 0, i.e. 
the original dilaton <fi is constant and we get the standard (flat) FRW - type cosmology in 
the radiation dominated era, a ~ t 2 /( Ar + 1 ). An important point is that this is an attractive 
solution, i.e. all solutions (with initial tp < 0) will eventually approach it ! So at late 
times we get the standard radiation dominated era of cosmology without any fine tuning. 
If a < —\/N the attractive solution will be f/l = —\fN. 
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Figure Captions 

Fig. 1. The solid curve represents the adiabatic variation of energy E as a function of 
scale A. Note the duality symmetry A — > —A. Near the Planck scale (A ~ 0) the 
energy is more or less independent of scale, as we are in the Hagedorn regime 
of string thermodynamics. For large radii, we enter the radiation dominated era 
and E oc exp — A. The dashed curve represents the effective energy if the winding 
modes (and momentum modes for A << 0) fail to annihilate. 

Fig. 2. The 2d black hole geometry is represented here. The cosmo logical solution cor- 
respond to regions II and IV. If we wish to deal with a compact universe, we 
are limited to the wedge drawn here. In this case the curve drawn in region iT 
inside the wedge represents the space, with its endpoints identified. 

Fig. 3. The 'accordion-like' Penrose diagram of 2d cosmology. An observer in universe n 
can end up in any universe m > n. The regions inside the curves represent the 
universe if we wish to consider compact space. 
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